Catalytic reactions with bulk-mediated excursions 
Mixing fails to restore chemical equilibrium 
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In this paper we analyze the effect of the bulk-mediated excursions (BME) of reactive species on 
the long-time behavior of the catalytic Langmuir-Hinshelwood-like A -\- B ^ Q reactions in systems 
in which a catalytic plane (CP) is in contact with liquid phase, containg concentrations of reactive 
particles. Such BME result from repeated particles desorption from the CP, subsequent diffusion 
in the liquid phase and eventual readsorption on the CP away from the intial detachment point. 
This process which leads to an effective superdiffusive transport along the CP. We consider both 
"batch" reactions, in which all particles of reactive species were initially adsorbed onto the CP, and 
reactions followed by a steady inflow of particles onto the CP. We show that for "batch" reactions 
the BME provide an effective mixing channel and here the mean-field-type behavior emerges. On 
contrary, for reaction followed by a steady inflow of particles, we observe essential departures from 
the mean-field behavior and find that the mixing effect of the BME is insufficient to restore chemical 
equilibrium. We show that a steady-state is established as f — > oo, in which the limiting value of 
the mean coverages of the CP depends on the particles' diffusion coefficient in the bulk liquid phase 
and the spatial distributions of adsorbed particles are strongly correlated. Moreover, we show that 
the relaxation to such a steady-state is described by a power-law function of time, in contrast to the 
exponential time-dependence describing the approach to equilibrium in perfectly stirred systems. 

PACS numbers: 05.50.+q; 64.60.Cn; 68.43. De ; 82.65. -|-r 



I. INTRODUCTION 

Catalytically-activated reactions play an important 
role in various processes in chemistry, physics and biol- 
ogy. Such reactions are involved, as well, in many indus- 
trial and technological processes, in which the design of 
desired chemicals requires the binding of chemically in- 
active molecules, which recombine only when some third 
substance - the catalytic substrate - is present 0, 0, . 

One of the simplest examples of such catalytically- 
activated reactions, which will be discussed here, is pro- 
vided by the so-called Langmuir-Hinshelwood scheme 
Qi 1^ 13 ' This reaction scheme involves two types of re- 
active species - an ^ and a B, which are spread in a 
gaseous phase in contact with a solid surface - a catalyst, 
may adsorb onto the surface at specific adsorption sites 
(at constant rates Q')^'^^), desorb from them back to the 
gas phase (at constant rate Q'^^'^^), and enter into the 
reaction 



A + B 



(1) 



at a finite reaction rate K, as soon as any two of un- 
like adsorbed species appear at neighboring adsorption 
sites. The reaction product P desorbs from the surface 
instantaneously and leaves the system. 

Within the conventional mean-field approach 0, 0) 13 j 
(in which one discards the correlations in particles' distri- 
butions on the catalytic surface), one gets, in particular. 



(A) 

in the simple limit Q^^^ = Q 
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and at low particles densities in the gaseous phase, the 
following large-i asymptotical behavior 



CA{t)=CB{t) 



Qa. 



ds 



K 
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(2) 



where CA,B{t) denote mean surface coverages by the A 
and B species at time t, respectively, while T deter- 
mines the characteristic time at which the value Coo = 
\/ Qads/K is approached. The expression in Eq.© can be 

readily generalized for arbitrary values of Q^^'^\ Q[^i^^ 
and for arbitrary particles densities in the gas phase, 
which will result in a somewhat more complex expres- 
sions for Coo and for the characteristic relaxation time T . 
Note, however, that the long-time approach to Coo will 
be still described by an exponential function of time. It 
is also important to emphasize that the state approached 
as t — > oo is believed to be a true chemical equilibrium 
state, in which the A and B particles distributions on 
the surface are (5-correlated and Coo is independent of 
the kinetic parameters. 

Within the last two decades a considerable progress 
have been made in the theoretical analysis of the kinet- 
ics of non- catalytic reactions 0, H, IE S H IS Oi 
II1|i1I3|1^JJ«JJ*I^J1- Here, a remarkable phe- 
nomenon of stochastic segregation has been discovered 



Siiaiiiia III El Es 111111111111111113, 

and the effects of correlations and fluctuations in parti- 
cles spatial distributions on the reaction course have been 
elucidated, which are in a striking contrast with the con- 



2 



ventional mean-field picture '2^. These studies resulted 
in the inception of a novel interdisciplinary domain on the 
boundary between conventional chemistry and physics - 
the fluctuation-dominated chemical kinetics. 

Following the early works on the fluctuation phenom- 
ena in chemical reactions 0: S S 13 1 ZifF and collab- 
orators 0, 1^ have questioned the predictions of the 
mean- field approach in Eq. ^ for the catalytic Langmuir- 
Hinshelwood scheme. Focussing on the specific example 
of the oxidation process of the carbon monoxide on plat- 
inum surfaces, CO -I- O2 ^ CO2 + O, Ziff et al have 
observed a behavior which is by far richer and goes far be- 
yond than the traditional mean-field predictions. In par- 
ticular, they have discovered that as the CO gas pressure 
is lowered the system undergoes a first-order transition 
from a CO saturated inactive phase (zero rate of CO2 
production) into a reactive steady state (non-zero rate of 
CO2 production) followed by a continuous transition into 
an 02-saturated inactive phase. This continuous transi- 
tion was shown to belong to the same universality class 
as the directed percolation and the Reggeon field theory 
[23| . Different aspects of the kinetic and equilibrium be- 
havior in this model have been scrutinized, revea ling the 
importance of many-particle effects. [2^ l2^ 123 . |27l l2^ 

El El HI [34, 35, silalSEin^flMtEar 

An essential ingredient of the Langmuir-Hinshelwood 
model is that the phase confronting the catalytic surface 
and acting as a reservoir of particles is gaseous. This is 
not the case, however, in many instances. For many im- 
portant applications, especially in biological and chemical 
systems, the catalytic substrate appears to be in contact 
with a liquid phase which comprises concentrations of 
reactive species (see, Fig.l). 

In such systems the reaction kinetics and equilibrium 
properties may be affected by yet another important pro- 
cess, not included in the previous models ; namely, here 
the particles can performs long-range concerted excur- 
sions inside the bulk liquid phase. That is, , as depicted 
in Fig.l, the particles, adsorbed onto the surface being 
in contact with the bulk liquid phase, can desorb, dif- 
fuse rapidly (with the diffusion coefficient being several 
orders of magnitude larger than that for a surface dif- 
fusion) within the bulk phase, and then return to the 
surface at a new position far away from the detachment 
point. 

Indeed, it has been shown both experimentally and 
theoretically (see, e.g., Refs- liS li llisllM liTll and refer- 
ences therein) that such bulk mediated excursions (here- 
after abbreviated as the BME) can be the principal kind 
of motion for particles. In biological systems, such bulk- 
mediated excursions might play a significant role since 
here most of the systems include membranes surrounded 
by fluid environments [4^ ; as an example, one may con- 
sider the receptor-ligand reactions which take place on 
a membrane surface Is^l , or the catalytic efficiency 
of proteins, which cut the DNA molecules (one dimen- 
sional substrates) at specific sequences More gener- 
ally, BME may be involved in certain "searching" pro- 



cesses |52f . 

Several prominent features distinguish the situations 
with a liquid and with a gas phase in contact with a 
solid substrate. First, in the case when the solid surface 
is in contact with a liquid, the desorption of the adsorbed 
particles is generally much more pronounced than in the 
latter case ; hence, one expects that the BME process 
will be more frequent here. Second, appearing in the liq- 
uid phase after the desorption event, a desorbed particle 
will move diffusively, being multiply scattered by the sol- 
vent molecules. In view of the geometry of the system, 
here the motion relative to the surface is effectively one- 
dimensional, such that, after desorption and excursions 
in the bulk, any particle will be certain to return back 
to the surface, in contrast with the situation with the 
gaseous phase in which the particle may travel away from 
the surface almost indefinitely. 
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Solid Surface 

Fig. 1 - a sketch of the reactive 
Langmuir-Hinshelwood-like system with bulk-mediated 
excursions. Black and white circles denote A and B particles, 
respectively, (a) denotes a particle configuration, in which 
reaction takes place. Solid and dashed curved lines depict 
effective trajectories of A and B particles, respectively. 



Last but not least, repeated many times, the adsorp- 
tion/desorption events separated by the bulk-mediated 
excursions, will result in effective motion of any given 
particle along the surface. A most striking point here is 
that this motion is superdiffusive, such that with respect 
to its surface displacements, any given particle performs 
a "Levy walk" (see Refs.|54l ISa] for ample physical dis- 
cussion). Consequently, instead of a familiar Gaussian 
propagator, one finds [13, H EE EE El that here the 
distribution P{r, t) of particles' displacements along the 
surface is that of a two-dimensional Cauchy process and 
is characterized by a long 1 /r^ tail [sE] : 

Pir,t)^:^- r^^-^, (3) 
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where c = D/h, D being the particles' diffusion coef- 
ficient in the bulk liquid, while h — bQadslQ stands for 
the " adsorption depth" and b is the " capture range" - the 
distance over which a particle can directly be adsorbed 
in a single displacement step 0, [i^ • This 
implies, in turn, that due to the BME, the distribution 
6{r) of particles' displacements r along the surface obeys 
lkd4 i45kdS» 47] : 

4>{r)^^, r>r* = VDF, (4) 

t* being the typical time between the readsorption events, 
t* = D/{Qadsb)^ EIjSI- The impact of such a peculiar 
transport on the kinetics and equilibrium properties of 
the reaction process in Eq.Q has not been elucidated 
theoretically up to the present time. 

In this paper we analyze the effect of the bulk-mediated 
excursions of the reactive species on the long-time behav- 
ior of the catalytic Langmuir-Hinshelwood-like reactions, 
Eq. Q , in terms of a simplified model, the model captures 
only some basic features of the physical system, but still 
allows to draw several important conclusions. First of 
all, we concentrate here on a totally symmetric situa- 
tion, in which the mean concentrations of the reactive 
particles in the liquid phase, as well as their adsorption 
and desorption rates, are equal to each other. Further 
on, in this model, we focus on the events taking place 
on the catalytic surface and introduce the presence of a 
semi-infinite liquid phase containing concentrations of re- 
active species in an indirect fashion. That is, we assume 
that continuous inflow of A and B particles, dispersed 
in the semi-infinite liquid phase, to the catalytic surface 
can be modelled as a source, which creates A and B di- 
rectly on the surface, independently of each other and at a 
constant production rate Qads- Furthermore, we suppose 
that the BME can be taken into account by letting the 
adsorbed particles to perform random, long-range hop- 
ping motion along the adsorption sites of the catalytic 
surface with suitably chosen hopping probabilities, deter- 
mined by Eq.Q. Consequently, we propose here a two- 
dimensional model which includes two types of reactive 
species, which react upon encounters, perform long-range 
(Levy or, more spcifically, GiUis- Weiss H^E^) walks on 
the lattice and are continuously introduced onto the lat- 
tice from a reservoir maintained at a constant chemical 
potential. 

We hasten to remark that this model serves only as 
a first approximation of the real physical system and 
there are several other important processes, which may 
influence the kinetic behavior of the reaction process in 
Eq.Q. First of all, an assumption that the "intensity" 
of particles creation on the lattice is not varying in time 
may be inadequate. As a matter of fact, here, in view 
of the effectively one-dimensional geometry and diffusive 
transport in the liquid phase, non-homogenous particle 
density profiles in the direction perpendicular to the cat- 
alytic plane, characterized by a "depletion" zone, may 
emerge. This will result in the "intensity" decaying in 



time. On the other hand, here we overestimate mixing 
effect of long-range BME, supposing that for any parti- 
cle a jump on distance r along the lattice, once chosen 
with the probability distribution in Eq.Q), is executed 
instantaneously, while in reality the transport via BME 
on this distance takes some time interval, which is a ran- 
dom variable having a broad distribution. Consequent 
analysis of these effects requires much more complex ap- 
proach, which is currently being carried out |60| . 

Finally, we would like to note that, apart from its 
relevance to the reaction process in Eq.(^ for a cat- 
alytic surface in contact with a liquid phase, our anal- 
ysis sheds the light on the conceptually important ques- 
tion of the effect of mixing on the fluctuation-induced 
kinetics. As we have already remarked, for batch A + 
B reactions, i.e. reactions in which the particles 
of the reactive species are all initially introduced into 
the system, in case of equal mean particles densities 
and diffusive transport, diffusion appears to be a non- 
efficient mixing channel and the like species tend to 
segregate spontaneously in the reaction course, which 
causes deviations from the text-book kinetic behavior 

iSiliailiailllillillillillilliliaiil. In 

presence of a steady inflow of reactive species, this effect 
gets dramatically increased and the deviations from the 
conventionally ex pec ted behavior |23 are getting even 
more pronounced ]TlL ITllTlIll 11611171 l39l| . On the other 
hand, it has been shown recently in Refs. |6lL Ig^ . Is^ . 
which analyzed kinetics of the batch A + B ^ reac- 
tions, involving particles which execute Levy walks, that 
in this case, under certain conditions, accelerated diffu- 
sion destroys effectively the particle-particle segregation 
and the mean-field behavior prevails. One may now pose 
quite a legitimate question whether in situations with a 
steady inflow of reactive species the Levy walks would 
provide a sufficiently fast mixing channel and overcome 
the strong tendency for segregation of like species in the 
course of the process in Eq. ? Our answer is negative. 
We show that in the situation under study, despite the 
fact that we strongly enhanced mixing supposing that 
the BME are executed instantaneously, the state reached 
by the process in Eq. as i — > cx) is not a true chemical 
equilibrium but only a steady-state. As a matter of fact, 
we proceed to show that the long-time concentration Coo 
appear to depend on the kinetic parameters, such as, e.g., 
particles diffusivities in the bulk, particles distributions 
on the lattice are very strongly (algebraically) correlated 
and moreover, the long-time approach to such a steady- 
state is essentially delayed, as compared to the exponen- 
tial dependence in Eq. - it is described by a power-law 
function of time ! 

The paper is structured as follows. In section 2 we 
introduce the model and basic notations. In section 3, 
focussing on the case of batch reactions, we present our 
analytical approach and reproduce several known results. 
In section 4 we analyze the steady-state behavior in mod- 
els with steady particles input, which mimics Langmuir- 
Hinshelwood scheme with bulk- mediated excursions, and 
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discuss the long-time approach to such a steady-state. Fi- 
nally, we conclude in section 5 with a brief summary of 
our results and discussion. 

II. MODEL 

Consider a two-dimensional regular lattice which is 
brought in contact with a reservoir of particles of two 
types - an A and a B, maintained at constant chemical 
potentials and /is. Here we restrict our analysis to 
the special case ij.a = fJ-B- The particles of both species 
may adsorb onto the lattice at constant rate Qads , desorb 
from the lattice at rate Q, an event followed by a long- 
range instantaneous jump of distance r with probability 
(j){r) and an immediate re-adsorption. The particles then 
react at a constant rate K according to the scheme in 
Eq. Q as soon as two of them of unlike species appear on 
the same lattice site. In most of our analysis we will focus 
on the limit K oo, which will allow us to emphasize 
the "statistical physics", rather than purely "chemical" 
effects. We will discard the hard-core exclusion between 
like and unlike species, assuming that the particles' cov- 
erages are sufficiently small. 

Now, the long-range jumps performed by the parti- 
cles of reactive species will be described here within the 
framework of the Gillis- Weiss random walks (also re- 
ferred to sometimes as the Riemann walks l59l'l . which 
represent the lattice version of Levy flights v4\ in the 



limit |r| ^ 1. We will use here a bit more general defini- 
tion of 0(r), than that in Eq.l^J, and supporse that (j>{r) 
is given by 

= (5) 

Note that the distribution in Eq.© reduces to the one 
in Eq.|01) in the particular case when d = 2 and /z = 1. In 
this case, the parameter £^ — r* . Note also that with this 
definition of the elementary jump probability, the mean 
square displacement per step is infinite for all < 2, 
which implies that such a random walk has an infinite 
variance |l64i65j . The long-tailed distribution of the jump 
probabilities permits long range jumps and generates a 
super-diffusive regime. Gillis- Weiss walks lead to anoma- 
lous diffusion, associated to the dynamic exponent 2/fi 

for /i < 2, |r| ~ i^/'', and to conventional diffusion for 
/i > 2, corresponding to Gaussian random walks, ^ t. 
In the case of interest here, i.e. for fJ, = I and d — 2, 
which case mimics the reaction in Eq.(Q mediated by 
bulk excursions, one has that in regard to surface dis- 
placements, the particles execute random ballistic-type 
(with an infinite vellocity) motion with |r| ^ t. 

Let now C^(r,t) and CB(r, denote the local, (at 
point with vector r) , time-dependent coverages of A and 
B particles, respectively. Evolution of these properties is 
governed by the following rate equations 



CA{r) = -KCa{v)Cb{v) - - ^ 0(r' - v)Ca{v) + -Y^cj^iy - v')Ca{v') + Qi^](r, t) 



(6) 



Cb{v) = -KCa{v)Cb{v) ~-Y. - r)CB(r) + 1 ^ 0(r - v')Cb{v') + Qi^](r, t), 

Td , Td ^ 



(7) 



where the dot denotes the time-derivative, the first term 
on the rhs describes the decrease in particles' coverages 
due to the reaction events, the second and the third terms 
describe departures and arrivals of the particles at the 
site r at time t due to long-range jumps, respectively. 
Note that the summation in the second and the third 
terms on the rhs of Eqs.JHI) and (jZj) extend over all lattice 
sites, which signifies the long-range character of particles' 
migration, while denotes the time each particle typi- 
cally spends on each lattice site between the desorption 
events, Td — Q^^- In turn, the fourth terms on the rhs of 
Eqs.jnj and {Tj) describe the (random) contributions to 
particles' coverages due to adsorption of particles from 
the reservoir, which mimics, in our model, the presence 
of particles in the bulk liquid phase. In the usual fash- 
ion, we admit the following statistical properties of these 



"source" terms : 

(d(r,t))-(Qi^](r,i))-Qa.., (8) 

+ A, i + r)) = Q^. + Qads5{X)5{T), 

(9) 

{q[2 (r, + A, i + r)) = Q^. + Qad.<5(A)<5(r), 

(10) 

and 

(Qi^l(r,OQlfi(r + A,t + r))=0, (11) 

where ^(A) is the delta-function, A is the correlation pa- 
rameter and the angle brackets (. . .) , here and henceforth. 
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denote the volume averages. Consequently, we stipulate 
that creation of particles on the catalytic surface proceeds 
completely at random, in space and in time, at a con- 
stant rate, which describe arrival of the particles located 
initially in the bulk liquid phase at progressively longer 
distances in the direction perpendicular to the catalytic 
surface. 

In what follows, we discuss the behavior of the solu- 
tions of the dynamic rate equations, Eqs.(|HJ and un- 
der different physical conditions. 

III. BATCH REACTIONS 

To set up the scene, we discuss first the effect of the 
BME on the kinetics of A + B in the "batch" 
reaction case ; namely, in situations in which all par- 
ticles of the reactive species which were dispersed in- 
tially in the bulk liquid phase were asorbed onto the sur- 
face by some "rapid quench" ; steady inflow of reactants 
by the external source is supposed to be absent here, 

'3ads('"'*) = Qads = ^- supposc, as wcll, that at 
t > particles' desorption from the catalytic surface, and 
consequently, the BME, are allowed. 

In this case, we assume that the initial particles' distri- 
butions on the lattice are random Gaussian, J-correlated 
with mean coverages Cq ; that is, CA(r,0) and CB(r,0) 
obey : 

(C^(r,0)) = (CB(r,0)) =Co, (12) 
(^CA{r,0)CA{r + \,0)) =C^ + Co5{X), (13) 
(^CB{r,0)CB{r + X,0)) ^C^+CoS{X), (14) 

and 

(c^(r,0)Cs(r-l-A,0)) =0. (15) 

Now, to analyze the time evolution of the mean par- 
ticles' coverages, we make use of the analytical ap- 
proach first proposed in Ref.|^ for the description of 
the fluctuation-induced kinetics of irreversible diffusion- 
limited A + B ^ reactions. In this approach the hier- 
archy of the reaction-diffusion equations for the higher- 
order correlation functions has been truncated at the 
level of third-order correlations. Subsequent works (see 
Refs.milllll) generalized the approach to more com- 
plex reaction schemes, e.g. to reversible reactions or re- 
actions involving interacting particles, and also showed 
that such a truncation is tantamount to the assumption 
that fields CA,_B(r,i) remain Gaussian at all times; this 
implies that the fourth-order correlations decouple au- 
tomatically into the product of the pairwise correlations 
which insures that the third-order correlations vanish. 



Following Ref.|5'|, we write first the local coverages 
Cyi,_B(r,t) in the form : 

where 5CA,B{^,t) denote local deviations from the mean 

coverages C{t). By definition, ^(5CA,B(r, t)^ = 0. 

Further on, we introduce the pair-correlation func- 
tions : 

GAB{X,t) = (^6CA{r,t)6CB{r + X,t)y (17) 
GAA{X,t) = (^6CA{r,t)6CA{r + X,t)), (18) 

and 

GBB{X,t) = (^SCBir,t)SCB{r + X,t)y (19) 

A being the correlation parameter. Note that since we 
have assumed a totally symmetric situation with regard 
to the adsorption/desorption rates, the correlation fuc- 
ntions GAA{X,t) and GBB{.X,t) are obviously equal to 
each other at any time moment, G^^(A,t) ~ GBB{X,t). 
Next, averaging Eqs.© and {7|), we obtain 

C{t) = -K[C^{t)+GAB{t)], (20) 

where we have used the notation G'ab(0 — GABiO,t), 
i.e. Gab(0 is the value of the particle-particle pairwise 
correlations at distance A which is equal to the reaction 
radius R (here, in the lattice version of the model, we 
assumed that reaction takes place when two particles of 
unlike species appear at the same site, i.e. R = 0). 

Note now that Eq. (|20|l shows that the time-evolution of 
the mean coverage is ostensibly coupled to the evolution 
of the pairwise correlation function. Note also that if the 
correlations are supposed to be insignificant, GAB{t) = 
0, as one generally takes for the mean-field approach, 
one obtains from Eq.lHOJ) that C{t) = -KC^{t), i.e. the 
conventional " law of mass action" . This law yields C (i) ~ 
1/ Kt for systems of any spatial dimension and regardless 
of the way how the particles move in the system. On 
the other hand, assuming perfect, instantaneous reaction 
with K = oo, one gets that 

C{t) = V-GAB{t), (21) 

which represents the mathematical formulation of the 
segregation effect ; as a matter of fact, Ea. (|21|l shows that 
the time-evolution of the observable - the mean cover- 
age, is guided at any time t by the time-evolution of the 
pairwise correlations in the system ! In case when both 
species move diffusively, Eq.ETl entails an unusual ki- 
netic law C{t) ~ l/f^^^ which thus pre- 
dicts for d < 4 a slower time-evolution of the mean cov- 
erage than that defined within the mean-field approach. 
Note also that for finite K equation H21|) holds for times 
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t sufficiently large such that the fluctuation-induced law 
C{t) ~ l/f^^^ determines the long-time asymptotic be- 
havior. 

Now, to determine analytically the time-evolution of 
the mean coverage in our case with BME, we have to 
evaluate the time dependence of Gab (t) which embodies 
all necessary information on the initial fluctuation spec- 



tra and particles dynamics. Making use of Eqs.® to H20|l . 
neglecting the third-order correlations and employing an 
evident symmetry condition between the AA and BB 
correlation functions, {GAAi^,t) = GBB(A,t)), we find 
that GAB{^,t) and G'^^(A,t) obey the following system 
of reaction/transport equations : 



J 



and 



Gab{\ t) = -2KG{t) \Gab{\ t) + Gaa{\ i)] + -Y] c^{y)CxGAB{\ t), 

y 



(22) 



GAAiX,t) = -2KG{t) GAB{\t)+GAA{\t) 



Y^c^{y)CxGAA{\t), 



(23) 



where 



CxGab{\ t) = Gab{\ + y,t) + Gab{\ -y,t)- 2Gab(A, t). 

I 



(24) 



Note that Eas. lf^ . lf^ together with Eq.jSDj) are closed and 
with respect to the mean coverages and pairwise correla- 
tions, and permits the evaluation of these quantities. 

In order to solve the system of Eas. (|22|l . (|23|) and 
Ea. (|20|l . it is expedient to introduce a pair of discrete 
Fourier transforms : 



F(k,t) ^^i^(r,t)e*(''"-) 



and 



where B denotes the first Brillouin zone. 
Transforming Eas. H22|l and H23|) . we get 



GAB(k,t) = ^2KC{t) GABO^,t) + GAA{Kt) 

2 

Td 



GAAiKt) = -2KC{t) GAi3(k,<) + GAA(k,t) 



-{m-m]GAA{Kt). (26) 

Td 



in which equations we have made use of the condition 
(/)(— k) — 0(k), since random jump process under consid- 
eration is symmetric invarient. Equations H25(l and H26() 
are accompanied by the initial conditions, which follow 
from the ones in Eas. ((T^ to ((T3|l . 



GAA(k,0) =GBB(k,0) = Go, GAB(k,0) = 0. (27) 



^ _ N _ Solution of Eqs.(|2SI) to l(?7jl can be readily obtained ex- 

-I- — (^0(k) - (/)(0)jGAB(k,t) (25) plicitly and reads : 



Gab (k, = - ^ exp 



2(0(0) -0(k))l 



r 



1 - cxp ( - 4ii: / C{t')dt' 



(28) 



We focus first on the perfect reaction case, when K — 00. 
Here Gab (k, t) is simply 



Gab{K t) 



Go 



■ exp 



i(0(O)-0(k))-l, (29) 



and GAB{t), which enters Eg. 1)21(1 and governs evolution 
of the mean coverage, is given by 



GABit) 



(2n) 



GABiKt)dk. 



(30) 
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It is well-known (see, e.g. Refs.|59|, |6l|) that the leading 
in the small-/c limit behavior of structure function (l)(k) = 



cxp(— ir • k) follows : 

m - ^(k) ^ A|kr 



(31) 



where A is a constant, A — Tr£^/T{d + fi) sin(7r/i/2), V{x) 
being the Gamma-function. Consequently, for d — 2 the 
long-time behavior of the correlation function GAsit) 
obeys : 



GABit) = - 



Co 



2(27r)2 
Co 
^2(2^) Jo 

r (^^ 
-a^iOo — 



exp 



exp 



2 At. 

Td 

2At, 

Td 



ikr 



dk 



kdk 



(32) 



where is a dimensionless constant, = 

r(2//^)r2/''(2 Ai)sin2/^(7r^/2)/4i+i/^'/i7ri+2A'. Hence, 
the long-time decay of the mean coverage C(t) is given 
by 



(33) 



Generalization of the result in Eg. 1)33(1 for arbitrary d is 
straightforward and yields 



(34) 



i.e. precisely the result obtained using a different theoret- 
ical approach in Refs.[^|H3j which studied the kinetics 
of batch A + S ^ reactions involving particles per- 
forming Levy walks. Note also that for /i = 2, which 
corresponds to the case of a standard random walk, we 
recover from Ea.lt^^the celebrated fluctuation- induced 

law Git) ~ t-'^i^ EElEllHIiD!. 

On the other hand, for 
sufficiently small values of such as /i < d/2, fomrula 
(ESI would give a decay faster than the usual mean field 
1 jt law. Such a behavior is of course specific of the per- 
fect reaction limit K = oo. Actually the system becomes 
perfectly stirred for such small /x, and the segregation 
phenomenon seizes to exist, so that the usual 1/i decay 
is recovered. Similar effect has been also obtained for re- 
actions in inhomogeneous systems in Ref.^3|. 

We turn now to the borderline case d — 2 and /i = 1, 
which is of special interest here, since it corresponds to 
the reactions mediated by bulk excursions. Here, Eq.JHSJ 
entails the behavior C{t) ~ t~^, i.e. the decay of the 
mean coverge mediated by the Cauchy walks proceeds 
exactly in the same fashion, as the decay obtained within 
the conventional mean- field approach 20 1. This circum- 
stance has prompted the authors of Refs.^l|,(63| to con- 
clude that that such a long-range transport may serve as 
an effective mixing channel which suppresses effectively 
the "undesired" segregation effect. In the next section, 
we will examine whether this conclusion remains valid in 



case of reactions followed by a steady inflow of reactive 
species. 

Finally, we note that this special case, in which the 
" mean- field" , purely chemical component, and the de- 
cay of correlations contribute to the overall kinetics at 
the same rate, can be viewed from a different perspec- 
tive. Following the celebrated CoUins-Kimball treatment 
of imperfect diffusion-limited reactions 67] , one may stip- 
ulate that here the time-evolution of the mean coverage 
obeys effectively the second-order rate equation of the 
form 



Git) = -KappC^it) 



(35) 



where Kapp is the apparent or effective rate constant, de- 
pendent both on K and on the parameters of the Cauchy 
process. 

To determine Kapp, we proceed as follows : We no- 
tice first that for finite K the mean coverage defined by 
Ea. (|20|) cannot decrease faster than the t~^ law expected 
from the mean- field kinetics. This implies that the inte- 
gral J^G{t')dt' is divergent for t oo, and hence, at 
sufhciently large times 



exp 



-4:K I G{t')dt' 
'o 



« 1, 



(36) 



and Ea. (|32|l describes correctly the long-time evolution 
of the pairwise correlation function for finite values of K. 
Now, from Ea. (|32|l we have that in this case 



GAB{t) 



(37) 



Substituting this expression into Eg. 1)20(1 and searching 
for the solution of the resulting equation in the form 
G{t) = 1/Kappt, we obtain 



app 



2Tr^e 
CorlK 




3(:2 



(38) 



Recollecting next that Td ~ Q ^ and ^ = r* = V Dt* , 
where t* = D / Q^^Jd^ , we may rewrite Kapp as 




^ ^ Cob2 f KQads I _^ 



7r3 V DQ 



_ 2^' ( DQ 

(39) 

which relates the kinetic behavior of the "batch" reac- 
tions with the BME to the physical parameters describ- 
ing our model. 



IV. REACTIONS WITH STEADY INFLOW OF 
SPECIES 

Let us now consider the case when the reaction process 
in Eq.JQ) is accompanied by a steady inflow of particles 
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onto the lattice, which mimics in our model of reactions 
with bulk-mediated excursions arrivals onto the surface 
of particles initially dispersed in the liquid phase from 
progressively larger and larger distances in the direction 
perpendicular to the surface. The statistical properties of 
external sources are defined in section 2, Eqs.JHJ to 

Averaging Eqs.@ and {Tj), we find that in this case 
the time evolution of the mean coverages is governed by 



C{t) = -K [C{tf +GAB{t)\ 



cads- 



(40) 



On the other hand, we find that the pairwise correlation 
function GAB(A,t) is still determined by Ea.l|^. while 
GyiA(A,t) obeys (see, e.g. Refs. 0,123 for the details of 
the derivation in the diffusion-controlled case) : 



J 



Gaa{\ t) = -2KC{t) \Gab{\ t) + Gaa{\ t)] + -V] ^{v)CxGaa{\ t) + Qads 

I 



(41) 



We note, parenthetically, that discarding correlations, i.e. 
setting in Ea. (|40|l formally GAsit) = 0, we recover the 
mean- field result in Eq.Q, which claims that : i) the 
state approached as i — > oo is a true chemical equilib- 
rium, ii) Coo = a/ Qads/K and iii) Coo is approached 
exponentially fast. On the other hand, there is an am- 
ple evidence that all three of these statements are wrong 
because of the segregation effects in case when reactive 
particles perform conventional diffusive motion on a two- 
dimensional surface in presence of a steady inflow 0,113 ■ 
However, as demonstrated in Refs. 61, 63] and in the sec- 
tion 3 of the present paper, for batch reactions with the 
BME, which ultimately result in random ballistic-like 



J 



motion along the surface, the segregation effect is sup- 
pressed and, apart of some renormalization of an the re- 
action constant kinetic behavior follows an essentially 
mean- field type dependence G{t) ~ 1/t. Consequently, it 
is not a priori clear whether GAsit) governs (or even 
contributes to) the long-time evolution of the system in 
the case d = 2 and ^ = 1. 

To answer this question, we turn to the time-evolution 
of GAB{\t), defined by Eqs.(j23l and (gTJ. Applying the 
Fourier transformation, we find, after elementary calcu- 
lations, that in the long-time limit (in which the asymp- 
totic behavior in Ea. (|31() holds) the Fourier-image of the 
pairwise correlation function G^B(A,t) reads : 



GABiKt) 



^ads'^d 

"4A|kp 



1 — exp 



2A|k|' 



Td 



r 



exp 



2A|k|^ 4K / G{t")dt" 



dt'. 



(42) 



Now, some analysis shows that due to the presence of the 
function AK J^, C{t")dt" in the exponential, the second 
term on the rhs of Ea. H42l) is negligibly small compared 
to the first one. Consequently, we focus on the behavior 
of the dominant contribution. 

Consider first the long-time evolution of Gab (t) , which 
enters the rhs of Ea. H40() and hence, may affect the evo- 
lution of the observable - the mean coverage. Inverting 
the Fourier transform, we find then that in the long-time 
limit GAsit) obeys, in d-dimensions and for arbitrary /i, 

exp 



GABit) = - 



Qads^d 

4A(2^ 



1 



2A|k|^ 



Td 



sTd 



27r 



dk 



J^fi — d-\-l 



1 — exp 



2A|k|' 



Td 



which signifies that in such situations Cyis(i) is a grow- 
ing with t function. This is precisely the case discussed 
in Refs. \vA l39l |. which concerned with the behavior of 
two-dimensional A + B reactions involving diffusive 
species in presence of external uncorrelated input of par- 
ticles. On the other hand, for fi < d the stationary value 
GAB{t = oo) is finite. 

Consequently, in the case of interest here, i.e. for d = 2 
and /i = 1, we find from Ea. H43|l the following behavior : 



(43) 



Analyzing the behavior of GAB{t), defined in Ea. (|^ . 
one notices first that GAsit = oo) = oo for fi > d, 



GAB{t) 



Qads'^d , Qads^d 



AA 



WttAH' 



(44) 



which predicts that the stationary value of GAB{t), 
i.e.—QadsTd/^A is approached as a power-law, at a rate 
which is proportional to the first inverse power of time ! 
Note also that GABif) < 0, which means that pairwise 
correlation slow down the forward reaction. 
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Substituting the expression in Eq. H44|l into Eq. (|40|) and 
solving the resulting equation in the limit t oo, we 
find eventually that at sufficiently long times the mean 
coverage obeys 



which mimics the presence of particle's concentrations in 
the bulk liquid phase, the bulk mediated excursions fail 
to establish chemical equilibrium. 



C{t) = 



K 2t:DQ 



(45) 



where we have made use of relations ^ = r* and = 
Q"^, and the characteristic relaxation time tchar is given 

by 



^char 



2TrDQ J 



(46) 



Now, several comments are in order. First of all, we no- 
tice that in contrast to the behavior in Eq.©, predicted 
by the mean-field approach, the actual approach to the 
t = oo-state is algebraic, proportional to the first inverse 
power of time. This means that the long-time relaxation 
of the mean coverage to its t — > oo value is governed en- 
tirely by the time-evolution of the pairwise correlations. 
Second, Coo is generally different of the mean-field value 
{— \/ Qads/ K) and reduces to it only when the particles 
diffusion coefficient in the bulk liquid phase D ^ oo. The 
fact that Coo depends on such a "kinetic" parameter as 
D is also quite a prominent feature - it shows unambigu- 
ously that in the reaction process under study there is no 
true equilibrium but rather a steady — state ! 

To emphasize this point, we turn finally to the analysis 
of the particle-particle correlations on the surface in the 
limit t = oo, embodied in the pairwise correlation func- 
tion GAs(r) = G^B(r,t = oo). From Ea. l|l^ we find 
then that these stationary correlations obey for = 2 
and ^ ~ 1 : 

GAB(r) = -%^^ / dO [ dk eyi-p[ikr cos{9)], 
Wn^A 7o Jo 

(47) 

which yields, in the limit |r| — > oo, the following behavior 



GAB{r) = - 



ds 



1 



Sn^QD |r| 



O 



,13/2 



(48) 



Therefore, the decay of correlations in the particles' dis- 
tributions on the surface is algebraic, i.e. the particles' 
distributions show a quasi-long-range order. As a conse- 
quence, despite the fact that the BME process effectively 
mixes the system in the case of batch reactions restoring 
the mean-field-type behavior and suppressing the segre- 
gation effects, in the case with a steady inflow of particles, 



CONCLUSION 



To conclude, in this paper we have analyzed the ef- 
fect of the bulk-mediated excursions of reactive species 
on the long-time behavior of the catalytic Langmuir- 
Hinshelwood-like reactions in systems in which a cat- 
alytic surface confronts a liquid phase, containg concen- 
trations of reactive particles. Such bulk mediated excur- 
sions result from particles desorption from the lattice, 
subsequent fast diffusion in the liquid phase and even- 
tual adsorption on the surface far away from the in- 
tial detachment point. Repeated many times, such BME 
yield an effectively long-range particles transport along 
the catalytic surface with superdiffusive related proper- 
ties. We have considered both "batch" reactions, i.e. re- 
actions in which all particles of reactive species were ini- 
tially adsorbed onto the surface, and reactions followed 
by a steady inflow of particles onto the catalytic surface. 
The latter situations, under certain assumptions, mimic 
the presence of particles concentrations in the bulk liq- 
uid phase which act as a reservoir of particles. We have 
shown that for "batch" reactions, in accord with previ- 
ous analysis 0, |H3 , the BME provide a very effective 
mixing channel which suppresses the segregation effects, 
such that the mean-field-type behavior prevails. On con- 
trary, for reactions followed by a steady inflow of parti- 
cles, we observe essential departures from the mean-field 
behavior and find that the mixing effect of the BME is in- 
sufficient to restore chemical equilibrium. We show that 
here a steady-state is established as t ^ oo, in which 
the limiting value of mean coverages of the catalytic sur- 
face depends on the particles' diffusion coefficient in the 
bulk liquid phase and the particles' distributions on the 
lattice are strongly, algebraically correlated. Moreover, 
the relaxation to such a steady-state is described by a 
power-law function of time, in contrast to the exponential 
time-dependence describing the approach to equilibrium 
in perfectly stirred systems. 
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